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COMPUTER SIMULATIONS OF INTERACTION EFFECTS 
IN SOLITON TRANSMISSION, FOR STUDENTS TRAINING 


Andrei D. PETRESCU!, Andreea Rodica STERIAN’, Paul E. STERIAN* 


Rezumat. Articolul de fafa isi propune sa prezinte unele rezultate numerice ale 
simularii propagarii tri-solitonilor, bazate pe ecuatia Korteweg-de Vries (KdV), 
utilizand (Maplel1), un program puternic, ce permite realizarea calculelor 
numerice, trasarea si animarea reprezentarilor grafice, operarea cu expresii 
analitice. Este discutata interactiunea tri-solitonilor in prezenta perturbatiilor, 
pentru a gasi modalitatile de crestere a capacitatii de transmisie si de micsorare a 
ratei erorilor la nivel de bit. Aceste simulari au fost gandite pentru a fi utile atat 
proiectantilor sistemelor de transmisii digitale de date si studentilor, pentru a sigura 
o mai buna intelegere a acestor fenomene. Acest articol extinde unele rezultate 
obtinute pentru solutile uni si bi-solitonice ale equatiei KdV. 


Abstract. This paper aims to present some numerical simulations of tri-soliton 
propagation, based on Kérteweg-de Vries (KdV) equation using (Maplel1), a 
powerful program that permits to perform numerical calculations, plot or animate 
functions and manage analytical expressions. We discuss the tri-soliton interaction 
in perturbations’ presence, in order to find the modalities for raising the 
transmitting capacity and decreasing the bit-error rate. These simulations are 
thought to be useful both for the designers working in digital data transmission and 
students, for better understanding of those phenomena. This article extends some 
results for the one and two-soliton solutions of KdV equation. 
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1 Introduction 


In [1] we will discuss a model for solving Optoelectronic problems by means of a 
computer environment (Maple10) which allows both numerical and symbolic 
solving of a wide range of applications (2D or 3D plots and plots animation). 


In this paper we refer to the solitonic solutions of the Korteweg-de Vries equation 
(KdV) and we show some numerical and graphical analysis of the solutions of 
the KdV equation. 


Finally, we present and we discuss some 3D graphics and 2D animated graphs for 
the tri-solitonic solution of KdV equation. 
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2 Solitons in Optical Fibers 


As we know, solution of the nonlinear evolution equation is the localized 
nonlinear wave called soliton ({2], [3]). The existence of such a solution is due to 
the balance between fiber dispersion and nonlinearity ((4], [5], [6]). This 


phenomenon occurs in ordinary glass optical fibers, for 1 > 1.30 um, when the 
intramodal dispersion is positive ([7]). 


40 


In these fibers, the characteristic distance L of soliton propagation is ([3], [8], [9]): 


100 


Fig. 1. The characteristic distance of soliton propagation versus the soliton’s power 


7p =200 Wm (1) 
P 9 


where P is the peak power of the soliton. So, for P = 50 mW, L = 10 km; 
obviously, we may increase L if we use smaller peak power, as seen in fig. 1. 


If EDFAs with dispersion-shifted fibers are used, the transmission distance may 
be 3,000 km at a 40 Gb/s bit rate ({9]). 


3. The Fundamental Soliton, as a Solution of the NLS Equation 


Wave propagation in an optical fiber with nonlinearity can be described by means 
of the nonlinear Schrédinger equation ([3], [5], [6]): 


WE 1 OE = m o (F; + Fp) 


lin 


(2) 
cor °° or? : 


where Proniin 18 the polarization due to the nonlinear Kerr effect ([7], [9]). We 
suppose E is a wave packet, described by the expression 
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E(z,t) = A(z.) x expl j(@ot —Byz)]. (3) 
After some calculations, we will get the following expression: 
(4) 


0A | 1 OA BOA 


=—jy|Al A. 
a va 20 ill 


Normalizing this expression, we get the canonical form of the NLS equation: 


2 (5) 
OU =i: crepe), 


GAD iN? 7? 
The fundamental soliton, described by the equation: 


es Zep i 


1 1 (oS iy, oe (6) 
exp! S| oa | | a eXp| ——_| f—-_— 
i, v, T, v, 
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(7) 


verifies the nonlinear Schrédinger equation, as we can see by direct calculations: 


because tanh(7’) (7') 
oU OU 2 oU LU 
— =-U/ anh = —2|| Ui and — = j—; 
ar (02 \ 2H Bi oD (8) 
oU 10°U 2 
23 j\ iv py 
then : 


Using Maple11, we can visualize the 3D envelope of this solution (Fig. 2). 


Fig. 2. The 3D plotting of the envelope (7). 
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4. Korteweg—de Vries Equation Model 


Let us consider the Korteweg-de Vries equation (KdV): 
u(x,t) a t) oe & u(x,t) 
Ot a 


his one dimensional solitonic solution ([4], [5], [12], [13]), obtained and verified 
by using Maple11 program: 


+ 6u(x, 1) ——— =0 (10) 


vie ul 
OU) = ee ey 
cosh[k(x—4k‘t)] 
In fact, using Maple11, we can write: 
3 


Kay = G me 1) +60, Ne u(x, d+ E = 
x 


u(x, 0 


2 
: 2p = | | cosh(k (=x + 4 k? ay 
"| cosh(k (-x +42 1)) 


+2e({2 ele 4 
OX | cosh( k (lpet 4 ath 


ie 2 a) | nsf BT [iter os ry) (13) 


3 
OX" | cosh(k (-x + 4K t)) 
so, the solution u(x, f) 1s: 


2h (14) 
wash Dr + 4 RB)’ 


We must observe that, from the solutions of the free particle Schrodinger equation 
and using Maple11 we can easily obtain the solutions of KdV equation as follows: 


pears (15) 
by means of the Wronskian formula 
beXvaln)aamale (16) 
where 
W=Wcosh(k, (x—4k, t)), sinh(k, (x-4k, 1))..y,) (17) 


is the Wronskian determinant composed of 
wis) (18) 
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and WS) (19) 
5, (20) 

means B =k, (x-4k, t) (21) 
for E<0 (22) 
and Eexiken(aeh 4k. 1) (23) 
0<E, (24) 


Maple11 gives us: 
Soliton := proe(w) local L; L := In(w); RETURN(simplify(2xdiff{ L, x $ 2))) end proc 
W=y(5) (25) 
where y(é)=cosh(é) Corresponds to a negative energy E of associated Schrodinger 
equation. E:=kh(x- 4k 1) (26) 


w :=cosh(k(x— 4k’ 1)) (27) 


Fig. 3. The 3D visualization of k dependence of the equation (11) 
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The Wronksi matrix is MW1 := [cosh(k (x—4 k? t))] (28) 
and its determinant reads: W1 :=cosh(k (x—4 k?1)) (29) 


5. One Soliton Propagation. Numerical Results 


Maple11 helps us visualize via 3D plotting the k dependence of the equation (11); 
for: 


a) k=0.5, te[—4, 4] and’ x €[-25, 25]; 

> plot3d(2*.25/ (cosh (.5* (x-4*.25*t)))*2, x = -25..25, t =-4..4); 
b) k=1, te[-1,1] and xe[-25, 25]; 

> plot3d(2/ (cosh ((x-4*t)))*2, x = -25..25, t = -1..1 ); 
c)k=0.1, tef[-10,10] and x €[—25, 25]; 

> plot3d(2*.01/ (cosh(.1* (x-.04*t)))*2, x=-25..25, t=-10..10); 
d)k=0.3, te[-10,10]and x <[-25, 25]; 

> plot3d(2*.09/ (cosh (.3* (x+4*.09*t)))*2, x=-25..25, t=-10..10); 


we got the representations seen in figure 3 (a—d) and, obviously, we can study the 
pattern for whatever values of k. 

The program Maplel1l helps us plotting solutions (we can use animation or 
individual frames), for different values of the wave number k, for different 
moments ¢ and for an arbitrary interval for x; we can study the dependence of 
speed and amplitude onto k. 

Let k= 0.6, te{-6,—2, 2,6} and x e<[-25, 25]; Maple11 gives (see fig. 4): 


Fig. 4. Soliton propagation for k= 0.6 and te {-6, =2,2, 6} : a, b, c, d respectively. 
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{—-6, — 2, 2, 6} 


and *€l-2>: 25]. Maplel1 gives (see fig. 5): 


Now, for k = 0.2, S 


Fig. 5. Soliton propagation for k= 0.2 and te {-6, 2, 6} : a, b, c, d respectively. 
As we can see, both soliton’s speed and amplitude are k dependent ([3], [12]). 


6. Two Soliton Propagation. Numerical Results 


For study the interaction between two solitons we need another solution of KdV 
equation ([1]). 
We can obtain and verify this solution of the KdV equation by means of Maple11: 


W= Wy, (5), WS, )) 


We consider in this case where: 


w,(&, ) = cosh(&,).».w,(&) = sinh(, ) 
Esk (x-4k, 1) & =k, (x-4k, 8) 
w, = cosh(k, (x-4k, 1) wy =sinh(k, (x- 4k, 0) 


2 2 
ee hk, (x-4k, t inh(k, (x-4k, 
The Wronski matrix: yy? -= Ceenela ees 2) 


sink, (x-4k, t))k, cosh(k, (x-4k, 1) k, 
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2 2 
. W2 :=cosh(k, (-x + 4k, t)) cosh(k,(-x+4k, t))k 
has the determinant: a : )) (h ( ; ky 
~ sinh(k, (-x+ 4k,’ 1) sink(k, (x +4k,'1)) &, 


and the corresponding two-soliton solution: 
2 2 
u2 = -2(k, cosh(k, (-x+4k, 1) k, —k, k, cosh(k, (-x+4k, t)) +k, ky 
4 2.2 4 Bt a 
+k," cosh(k, (—x +4 k, ty" 2h = cosh(k(@at4 kt) k,°) y ( 
=cosh(k, (-x 4k,” 1)) cosh(k, (-x+ 4k, 1))k, 


2 
+ sinh(k, (-x + 4k, 1)) sinh(k, (—¥+4k,’ ))k,) 


b 4 


-40 -20 20 40 -40 -20 20 40 


Fig. 6. The interaction of two solitons, for: k; = 0.5, kp = 1, 
t e{- 6, 2, 6, 10} 


We check with Maple11 that it is a solution of KdV equation: 


3 
& u(x, 0) + 6 u(x, f) (5: u(x, 0) + & u(x, | = 


The second solution is: 
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2 2 
u2 =2(-k, k, cosh(k, (x-4k, 1) +k, k, cosh(k, (x-4k, 1)) —k, k, 
4 2.2 4 Be a 
~ k," cosh(k, (x-4k, t)) +k, +cosh(k, (x-4k, 1) k) re 
cosh(k, (x- 4k,’ 1)) cosh(k, (x- 4k, 1) k, 


— sinh(k, (x- 4k, 1)) sinh(k, e497 DVR): (30) 


With these solutions, Maple11 helps us visualize the interaction between the two 
solitons (Fig. 6). 

We can see both solitons regaining their shape after collision (i.e. after they pass 
one through another). Obviously, is better to use the animation feature of the 
Maple11 program, modify the parameters and observe the changes. 


7. Three Soliton Propagation. Numerical Results 


For the tri-solitonic solution one may obtain and verify (by means of Maple11) 
the following solution: 


vst Va: aes) 


Jo 2:(b—c) 
i2 


(ea) we Ve:(x-2:te) J. vst Vv a:(x—2:ta) }- 
J2 v2 


S: tanh| YE =2: te) = tanh 2 (#=2: ta) : 
[yee mal LEE) - rea ama| Lee tal) 
(ab) [ont _- ‘t'b) J» +e oe fa) J 


[Ve tn Va (x—2: ta) | FF oo Vb (x—2: tb) \) 


v2 v2 
2:(a—b) 


a SE a SE 


2:(¢€—a) 2 


= Ve(x 2 re) | arta Let8—2 20) 


where a, b, c are real, positive constants. 
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Maplel1 permits us to make some useful plots (see fig. 7) — specially animated 
plots - choosing various values for the parameters a, b and c. 


t= 1.6667 


Fig. 7. Animated plots for the three solitons interaction. 
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8. Conclusions 


Those computer experiments are very suitable for proofing to students in optical 
engineering and designers of optical circuits the main features of soliton 
propagation in optical fibers, as follows, reproducing those from [1]. 


The fundamental soliton appears as solutions of NLS equation. 


Solitonic solutions may be obtained solving KdV equation. Solving equation by 
means of Maplel1 program is useful because is quick, but the solution must be 
verified by hand. 


The existence of solitons is due to the balance between the dispersion and the 
nonlinearity of the medium. 


Both soliton’s speed and amplitude are k dependent; this dependence may be 
easily observed by modifying k in equation (14) and (15). 


For designers working in digital data transmission, the simulation of the solitons 
propagation along optical fibers is useful and may be performed successfully with 
Maple11. 


The animated plots and the 3D plots also are some useful features of this 
computer program; especially: the collision of two or three solitons and the k 
dependence of a soliton propagation can be studied by this means respectively. 


Solitons scatter elastically; after the collision, solitons regain their original shape 
and velocity ([8], [9], [15]). 


The only remaining effect of the scattering is a phase shift (i.e. a change in the 
position they would have reached without interaction). 


For solitons keep invariant shape and size is in accord with the conservation laws, 
one can infer that a profound link exists between integrable models and the theory 
of solitons. 
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